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We report the first observation of bound-state proximity 
resonances in coupled dielectric resonators. The proximity 
resonances arise from the combined action of symmetry and 
dissipation. We argue that the large ratio between the widths 
is a distinctive signature of the multidimensional nature of 
the system. Our experiments shed light on the properties of 
2D tunneling in the presence of a dissipative environment. 



Tunneling and dissipation are ubiquitous phenomena 
in physics. A detailed understanding of their combined 
action would be highly desirable given the relevance 
of the problem for atomic physics, condensed matter 
physics, chemistry and biology However, the incorpo- 
ration of dissipative effects is by no means trivial. Due to 
limitations of the available analytical and computational 
methods, up-to-date descriptions are still restricted to a 
few manageable cases, the prototype situation involving 
a bistable potential in ID 

In this Letter we report the observation of novel as- 
pects of tunneling in 2D potentials and its interplay with 
classical dissipation. In experiments utilizing microwave 
dielectric resonators, we find that symmetry not only 
plays a crucial role while shaping the eigenstates of the 
system, but also influences the way they couple to the 
external environment acquiring a finite width. In the 
observed resonance multiplets, we find that one of the 
members is extremely sharp due to the symmetry of the 
configuration. The large ratios of the observed widths 
appear to be a peculiar consequence of the multidimen- 
sional nature of the system. 

The experiments were carried out using MgTi dielec- 
tric cylinders placed between two parallel copper plates, 
30 cm square, separated by a gap I = 6.38 mm (Fig. 1). 
The disks had diameter D = 12.65 mm and dielectric 
constant Sr — 16. After establishing input/output cou- 
pling to the near field of the resonators by inserting coax 
lines terminated by loops, measurements of the trans- 
mission amplitude as a function of the frequency were 
performed using an IIP8510B network analyzer. 

The eigenvalue problem of a single dielectric resonator 
can be solved analytically by regarding the system as a 
waveguide along the direction z orthogonal to the plates 
The entire field configuration can be derived from 



the knowledge of the longitudinal components {Hz,Ez} 
alone, that separately obey the Helmoltz equation: 



{^^ + e){H,{v),E,{v)} = Q. 



(1) 
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Here, r — [cf), p, z) in cylindrical coordinates and k = 
(w / c) denotes the medium wave number for a mode 
at frequency w = 27r/ (e^. = 1 outside the dielectric) . For 
perfectly conducting walls, boundary conditions require 
kz — pi^/l, p integer. We have verified through explicit 
measurements of the field profile that p > 1, and all 
modes are evanescent with a decay constant close to the 
expected value = i/fc^ — uj^/c^ Q|. A generic mode 
of the dielectric is classified according to its azimuthal, 
radial and vertical quantum numbers (to, n,p). If to = 0, 
the mode has cylindrical symmetry and can be either TE 
or TM. Hybrid HEM modes arise whenever to > 0. There 
is no TEM mode for a dielectric guide. The agreement 
beween the calculated resonances and the data is found 
to be within 2% for all the peaks 

The quality factors Q of the resonances are deter- 
mined by the observed widths 7 in the frequency do- 
main, Q = 7/7. For a single resonator, estimates of Q 
are possible by calculating the ratio between the energy 
stored per unit time and the average power dissipated 
Q. Owing to the localized nature of the field eigen- 
modes, losses due to the open boundary conditions at the 
edge of the plates are irrelevant, power dissipation being 
introduced by dielectric and conductor losses. Detailed 
calculations, which yield results consistent with the mea- 
surements, indicate that finite absorption in the metallic 
plates outweighs dielectric losses by at least a factor 5, 
thus providing the leading dissipation mechanism |^ . In 
an equivalent time-domain picture, this implies that the 
metal acts as an environmental decay channel for the 
bound electromagnetic modes, the coupling between the 
dielectric and the metal being proportional to the copper 
surface resistence. 

When two dielectric resonators are placed in proxim- 
ity to each other, each resonance splits into two. The 
doublets have the structure of a broad resonance at a 
lower frequency fi accompanied by a narrow resonance 
at higher frequency fh- This effect is particularly pro- 
nounced for TM modes {Hz = 0). The most notewor- 
thy example is the TMqh single-disk resonance found at 
9.45 GHz with Qq « 70, which splits into two peaks with 
Qh ~ 2400 and Qi « 50 for an edge distance d — 1.0 mm 
(Fig. 2). The narrow peak can be experimentally as- 
signed to an antisymmetric _E^-field configuration by es- 
tablishing electric field coupling with the pick-up antenna 
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and by probing the behavior at the mirror symmetry 
plane Q . The doublet splitting as a function of the disk 
separation Af{d) ^ fh — fi is displayed in Fig. 3(a). The 
splitting vanishes exponentially with d until the limit of 
noninteracting resonators is approached. The measured 
decay constant is in good agreement with the single-disk 
value Kr — 0.45 mm"-'^, as expected on the basis of semi- 
classical estimates in a tunneling regime where Krd > 1 
l^. The widths 7;, 7/j and their ratio "fi/jh are plotted 
in Fig. 3(b) and in Fig. 4 (Dots) respectively. Again, 
single-disk behavior is recovered for sufficiently large d, 
where 7/, 7/1 — > 70- For small separations, the width 7/1 
is highly suppressed, leading to the high Q-values noted 
above. As indicated by Fig. 4, a maximum ratio ^i/jh 
~ 50 is seen at d = 0.73 mm. It is very remarkable that, 
thanks to the proximity effect, Q's in the range of 10'^ 
are achievable without resorting to closed-walls cavities. 

A quantitative account of the above results can be only 
achieved by numerically solving Eq. (1) with the appro- 
priate boundary conditions. Even if the problem is sim- 
plified since the z-dependence is separable, an accurate 
calculation of the lineshape factors Q requires the com- 
plete knowledge of the electric and magnetic field distri- 
bution within the cavity volume. By referring to [Q for 
more detail on the full electromagnetic analysis, our pri- 
mary goal here is to gain simple qualitative insights. Let 
us focus henceforth on the TMon configuration. In the 
two-disk system, the splitting into modes of well defined 
parity is easily understood as a consequence of the pertur- 
bation introduced by the resonator-resonator coupling. 
The latter is known to take contributions only from the 
unperturbed evanescent fields Ez{r) — cos{kzz) Ez{x,y) 
of each resonator Q . Experimentally, the antisymmetric 
mode is found to be able to store an extra amount of 
electromagnetic energy compared to the symmetric one 
0. This takes place through an amplification of field 
components (e.g., Ez,Hz), which do not contribute to 
power dissipation, whereby the higher observed Q. 

The fact that the stabilisation of the antisymmetric 
mode only manifests at small separations suggests to pic- 
ture the phenomenon in terms of a collective effect arising 
when two discrete states (the bound single-disk modes) 
are coupled to each other and, in addition, to a com- 
mon environment (the metal) that renders them unsta- 
ble. Similar effects are found in quantum physics, where 
they require an appropriate modification of the standard 
Weisskopf-Wigner decay theory ||^. The correspondence 
between electromagnetic (em) and quantum mechanical 
( Qm) systems is well established for stationary problems 
. In particular, taking into account the boundary con- 
ditions at the dielectric surface, the component Ez{x,y) 
in the waveguide plays the role of the wavefunction tp in 
a 2D quantum mechanical system, the dielectric medium 
corresponding to a square potential well at a Gxed en- 
ergy p3[ . Accordingly, the two-disk system maps into 
a 2D tunneling problem. In the presence of losses, the 



damping of the EM field amplitude is usually accounted 
through a complex frequency / — i"f/2 whose imaginary 
part provides the time decay rate Despite the fact 
that, due to the different time-dependent equations of 
motion, the correct mapping to complex energies of quan- 
tum unstable states is a nonlinear relation of the form 
(/ ~ *7/2)^EM ^ (e — «7/2)qm, a QM configuration which 
is stabilized against decay will still be mapped into an 
EM non-decaying state. 

A simple argument supporting the stability of the 
antisymmetric state goes as follows. Let |L) , \R) de- 
note degenerate ket states localized in the left, right 
well respectively. The two levels are coupled to each 
other by a tunneling perturbation of the form Ht = 
-T[\L) {R\ + \R) T > 0, and to a common contin- 
uum of states with a strength Wl, Wr. It is possible to 
show that the coupling to the environment mediates an 
extra interaction between the two discrete states, which 
strongly affects the decay properties of the combined sys- 
tem and thereby its spectral response If Wl ~ W^,, 
one predicts that the symmetric combination [\L) + 
corresponds to a Lorentzian resonance line at frequency 
/s, whose width is twice larger than the width of each 
single level coupled to a continuum of the same strength, 
while the antisymmetric combination [|L) — \R)], found 
at frequency /a, is completely stabilized. This behavior 
can be regarded as the counterpart of Heller's predictions 
for the proximity effect based on a point scatterer model 
I p^ . It is worth mentioning that the overall effect of the 
extra interaction indicated above is to dress the tunneling 
matrix clement T with an additional imaginary part. 

The actual situation, where the width of the antisym- 
metric mode is limited by the dielectric losses, would be 
more adequately modeled by invoking two distinct envi- 
ronments. To complement the previous analysis, we also 
explored a phenomenological description of dissipation in 
terms of an effective non-hermitian Hamiltonian . By 
introducing a drastic approximation, we only consider an 
effective ID potential projected along the horizontal sym- 
metry axis ||lo[| . Within the theory of multidimensional 
tunneling, this is supported by the fact that the inter- 
action is dominated by instanton orbits between the two 
centers Thus, we use a potential energy function of 
the form V pot{x) = Vo + iVi (for |a;| > d/2 + D), iV2 (for 
d/2 < \x\ < d/2 + D), Vb + iVi (for < d/2), all param- 
eters being real numbers. The imaginary terms iVi, iV2 
account for the losses outside and inside the double-well 
region respectively. We assume \Vi^2 \ /Vb,6 <C 1. 

The real part of Vpot{x) is depicted in Fig. 4 (inset). 
We allow for the possibility of a height barrier Vb ^ Vq 
to effectively include corrections arising from the 2D na- 
ture of the problem. The presence of extra-contributions 
to the tunneling interaction, which are lost in the ID 
model, is simulated by a more transparent barrier. The 
noninteracting limit corresponds to d ^ 00. For finite 
d, even and odd states are generated, with eigenenergies 
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Ep = sp — ijp/2, P = S, A. If K = Kr + iiii denotes 
the inter-well wave vector, each complex eigenvalue has a 
structure involving exponentials e"'^'''*, convoluted with 
oscillating functions of Kid whose details depend on the 
state |L), \R) The signature of tunneling shows up 
through the exponential dependence of the energy split- 
ting, Ae = eA — £s ~ e~ '^''^^ F {Kid) . The oscillatory terms 
in F are responsible for a "rippled" structure of the split- 
ting decay, which is apparent in the data (Fig. 3(a)). 

The transcendental equations determining Ea and Es 
have been solved numerically for different sets of param- 
eters with both Vh — and Vb < Vq and the results 
compared with the experimental ones [Q. The leading 
exponential decay of the energy separation and the asym- 
metric small-distance splitting of the widths are correctly 
predicted. However, we find a major difference between 
the two models in their capability to reproduce the ex- 
ceedingly large ratio between the widths. In the simula- 
tions with Vb = Vq, we were unable to reach ratios larger 
than 4, regardless of the values of the parameters Vi, V2, 
mainly affecting the absolute range of the widths. This 
order of magnitude is in agreement with independent re- 
sults on symmetry splittings of resonances due to semi- 
classical creeping orbits In the reduced-height con- 
figuration, the ratio "fs/ja can be controlled over a broad 
range (up to 50) by varying Vh/Vb- Some representative 
behaviors are summarized in Fig. 4 for barrier opacity 
Vb/Vo — 1, 1/4, 1/9. Maximum stability of the anti- 
symmetric mode is reached at an intermediate distance 
d, corresponding to the ■js/ja -peak value. The better 
qualitative agreement attainable with the reduced-height 
model suggests that dimensionality effects also play a key 
role in the experiment. 

In order to confirm the conclusion that a judicious use 
of symmetry leads to a Q-sharpening, we carried out ex- 
periments on a 3-disk system, with the disks placed at 
the vertices of an equilateral triangle. Resonance triplets 
are observed. In analogy to the 2-disk system, we predict 
that modes tranforming antisymmetrically with respect 
to reflections in each of the vertical mirror planes show 
enhanced stability against dissipation Q . We focus on a 
single-resonator mode with m = 3 at 10.8 GHz, whose 
behavior is displayed in Fig. 5. A very sharp component 
at intermediate frequency is clearly seen. The Q factor is 
increased by roughly a factor 20 compared to the original 
one. The sharpening effect turns out to be very sensitive 
against symmetry-breaking effects. The influence of a ge- 
ometric symmetry-breaking has been studied by shifting 
one of the disks by b along the bisectrix of the triangle. 
It is evident from Fig. 5 (inset) that the sharp resonance 
is dramatically affected, the Q factor being exponentially 
degraded. No sensible change is observed for any of the 
broad resonances. 

In the spirit of the original definition by Heller JT2[ , 
our observations indicate that interesting proximity phe- 
nomena arise in the spectral response of nearby systems. 



Proximity resonances have been recently detected in the 
scattering of a TEM electromagnetic mode in a parallel- 
plate waveguide jl^]. Despite some superficial similarity 
with the present work, it is essential to realize that our 
experiments probed a completely different regime of the 
microwave field, where direct evanescent-wave coupling 
between bound modes rather than scattering resonances 
from two dielectrics illuminated by the same wave field 
were investigated. In particular, the power law behavior 
characterizing scattering states [|l2| should be contrasted 
with the exponential dependences that are intrinsically 
associated with tunneling. Thus, bound-state proximity 
resonances form a novel complementary manifestation of 
a similar physical phenomenon, whose detailed under- 
standing poses new challenges to both numerical simula- 
tions and semi-classical treatments. 

Our results have a variety of implications. First, the 
tunneling interaction in a 2D integrable potential has 
been probed sensitively. There is no difficulty, in princi- 
ple, to extend the experimental work to chaotic poten- 
tials where important results on complex periodic orbits 
theory have been obtained Second, the experiments 
demonstrate how symmetry properties can be usefully 
exploited to protect a system against the effects of its 
environment. This general mechanism provides a unify- 
ing explanation for proximity resonances, regardless of 
the unbound (scattering) or bound (confined) nature of 
the wave field. Third, our work can be related to recent 
observations of the symmetry splitting between optical 
modes in photonic molecules where, however, the 
behavior of widths was not addressed. It is conceivable 
that some counterpart of proximity phenomena may be 
relevant on the mesoscopic scale as well. From the prac- 
tical perspective, the possibility of symmetry-based Q- 
amplification in electromagnetic or photonic structures 
represents another exciting area of applications. Finally, 
the electromagnetic phenomenon evidenced here displays 
intriguing similarities with concepts investigated in the 
context of quantum dissipative processes jl^. The 
possibility of establishing some mapping between the 
electromagnetic and quantum realm even in the presence 
of dissipative mechanisms would clearly open up a fruitful 
arena of interchange and deserves further investigation. 
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analysis software and N. D. Whelan for discussions. 



[1] J. Jortner and B. Pullman, Tunneling (Redidel, Dor- 
drecht, 1986). 
[2] A. J. Leggett, Rev. Mod. Phys. 59, 1 (1987). 
[3] Dielectric Resonators, edited by D. Kajfez and P. Guillon 



3 




(Artech House Inc., Dedham, MA, 1986). 
[4] K. Pance and S. Sridhar, in preparation. 
[5] S. C. Creagh and N. D. Whelan, Phys. Rev. Lett. 77, 

4975 (1996). 

[6] S. M. Barnett and P. M. Radmore, Methods in Theoret- 
ical Quantum Optics (Clarendon Press, Oxford, 1997). 
[7] S. Sridhar, Phys. Rev. Lett. 67, 785 (1991). 
[8] L. Viola (unpublished). 

[9] S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge University Press, Cambridge, 1995). 

[10] Such a procedure amounts to replacing the physical circu- 
lar resonators with infinitely extended dielectric regions. 

[11] A. Wirzba and P. E. Rosenqvist, Phys. Rev. A 54, 2745 
(1996). 

[12] E. J. Heller, Phys. Rev. Lett. 77, 4122 (1996). 
[13] J. S. Hersch and E. J. Heller, Phys. Rev. Lett. 81, 3059 
(1998). 

[14] M. Bayer et al, Phys. Rev. Lett. 81, 2582 (1998). 
[15] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 
Lett. 81, 2594 (1998). 



FIG. 1. View of the experimental setup. Shown in the 
photo are two cylindrical dielectric resonators, and input and 
output antennas. 




FIG. 2. Experimental proximity resonances for two cou- 
pled resonators operating in the TMoii mode at 9.45 GHz. 




FIG. 3. (Top) Doublet frequency splitting A/ and (Bot- 
tom) Widths 7i,7h vs. distance d for the two-disc proximity 
resonance around 9.45 GHz. The result of an exponential fit 
in the region Krd > 1 is also shown. 
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FIG. 4. Ratio between symmetric and antisymmetric width 
vs. distance d/D for experimental data (Dots) and various 
implementations of the double-well potential shown in the 
inset. Units where h — 2m — 1 have been chosen. Vb = 900 
in units D-^ and Vi/Vo = -3 ■ 10-^ V2/V0 = -3 ■ 10"^ The 
barrier height is T4/Vo=l (a), 1/4 (b), 1/9 (c). 
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FIG. 5. Resonance triplet for three coupled resonators 
operating in the hybrid (3,1,1) mode at 10.8 GHz. (Inset) 
Quality factor of the sharp component as a function of the 
symmetry-breaking parameter b/D. 



